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This paper is devoted to robust, Predictor-based Model Reference Adaptive Control 
(PMRAC) design. The proposed adaptive system is compared with the now-classical Model 
Reference Adaptive Control (MRAC) architecture. Simulation examples are presented. 
Numerical evidence indicates that the proposed PMRAC tracking architecture has better 
than MRAC transient characteristics. 
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original system matrices 
extended system matrices 

Hurwitz matrices for reference model and predictor 

unknown diagonal matrix, with positive elements 

control input vector 

system state vectors 

controlled system output vector 

state and control dimensions 

matched uncertainty 

known regressor vector 

dimension of the regressor vector 

matrix of constant unknown parameters 

external bounded time-varying vector of commands 


I. Introduction 

T his paper is devoted to the design and analysis of a state-predictor based direct adaptive controller for Multi- 
Input-Multi-Output (MIMO) dynamical systems, with matched uncertainties. Our main intent is to evaluate 
possibility of improving MRAC transient characteristics during a tracking task. Towards that end, we propose the 
use of a state-predictor in formulating adaptive laws. Specifically, starting with a direct MRAC system, [1, 2, 4], we 
design a state predictor, while using the system full state measurements. Then through Lyapunov’ stability analysis, 
we propose a predictor-based modification to the MRAC laws. This modification constitutes the novelty of the 
results that are reported here. 


The rest of the paper is organized as follows. In Section II, MIMO system dynamics are defined and the control 
goal is formulated. In section III, predictor-based MRAC laws are derived, and stability proof is conducted using 
conventional (in adaptive control) Lyapunov-based arguments. Simulation example is presented in Section IV, 
followed by conclusions that are given in Section V. 


II. Model Definition 

In this section, we formulate the system dynamics, pose the control problem, and derive tracking error dynamics 
for adaptive control design. We begin by considering a class of MIMO uncertain systems in the form: 
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X, = A p x p + B P A « + ©rf (*„ ) (2-1) 

\ / 

where x p E:X p C is the system state vector, wEt/ G R m is the control input, 

d(x p )=e T d * d (x p ) (2.2) 

is the linear-in-parameters state-dependent matched uncertainty, Q d E R Nxm is the matrix of unknown constant 
parameters, and (x p )e 7?^ is the known N - dimensional regressor vector, whose components are piece-wise 
continuous functions of x„ . Also in (2.1), 7? E7?" Xffl is cons tant and known, A n E7?” x ” is constant and unknown, 
and A E R mxm is constant, diagonal and unknown matrix with positive diagonal elements. 

The control goal of interest is bounded tracking , that is one needs to define u such that the system controlled 
output 

y = C p x p EiR m (2.3) 

tracks any bounded and possibly time- varying command r(t)E:R m , with bounded errors, in the presence of the 
system uncertainties &’ A > e„}. 

Let 

e y( t )=y( t )- r ( t ) ( 2 * - 4 ) 

denote the system output tracking error. Augmenting (2.1) with the integrated output tracking error , 

e yl (0= f e y ( T V T ) [ e yi = — ) ( 2 - 5 ) 


yields extended open-loop system dynamics: 

x = Ax + B\(u + d ^ ))+# c t (2.6) 

where x = (e T yI x T p )e X C R n is the extended system state vector, whose dimension is n = n p +m . The system 
matrices (^4, B , Z? c ) are 


C P \ B 

®n p xm A p B p 


y = ( 0 pxm Cp)x = Cx (2.8) 

_ jy 

represents the system controlled output. T o summarize, in this paper we will consider MIMO uncertain systems in 
the form of (2.6): 

* (0 = A x (0 + B A (« + ® ,/ fa B c r (t) 

y{t)=Cx(t) 

with an unknown constant matrix AE:R nxn , known matrices BE:R nxm , B c EZ?" xm , CE:R mxn , an unknown constant 
diagonal positive-definite matrix A E R mxm , and an unknown matrix of constant parameters & d E R Nxm . The state- 
dependent regressor vector (x p )e R n is assumed to be known. The control goal is bounded tracking of any 

bounded time-varying reference signal r(t^E:R mxl by the system output y(t)E:R mxl , in the presence of constant 
parametric uncertainties (^4, A, ) . 


Assumption 2.1 (Model Matching Conditions) 
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Given a reference Hurwitz matrix A ref and an unknown positive-definite diagonal constant matrix A , there must 
exist a constant possibly unknown gain matrix K x such that 

A e/ =A + BAK T x (2.10) 

Using the assumed relation (2.10), we rewrite the system dynamics in the form: 

x = A, 4 x + B a(u - K T X x + @ T d (x p ))f5 c r (2.1 1) 

Control input u is chosen as: 

u = K T x x-& T d <S> d (x p ) (2.12) 


where K x (t)E:R nxm and & d (t)E:R Nxm are adaptive time-varying matrices, whose dynamics will be defined later. 
Substituting (2.12) into (2.1 1), yields: 

/ \ 


x = A ref x + B A 


(k x -K x )x-{® d -@ d )<S> d (x p ) 


+ Br 


= A^x + BA(AK T x x-\e T d <t> d (x p ))+B c r 
Based on (2.13), consider the following reference model. 

Kef = A refXref+B c r 

and let 


e = x - x, 


ref 


represent the system tracking error. Subtracting (2.14) from (2.13) results in the tracking error dynamics , 

e = + BA(AK T x x- (x p )) 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


Remark 2.1 

Based on (2.16) and using Lyapunov arguments, direct Model Reference Adaptive Control (MRAC) laws can be 
derived [1], [2] in the form 


K x =-r x xe T P^B 


(2.17) 


[e rf = r 0 , e T P refB 

where P ref = > 0 is the unique symmetric positive definite solution of the algebraic Lyapunov equation, 

A T re fP ref +P ef A ref =-Q ref (2.18) 

and Q ref =Q T ref >0. Relations (2.12) and (2.17) solve the tracking problem, with globally asymptotically stable 
closed-loop dynamics. This particular solution is valid for any symmetric positive definite rates of adaptation T x 
and r 0 ^ . However, if these matrices have large singular values, then unwanted transient oscillations will often 
occur in the system. 


III. Predictor-Based MRAC 

In this section, a PMRAC -based solution to servomechanism problem is presented, for the class of uncertain 
systems in the form of (2.9). Towards that end and based on (2.1 1), define the so-called predictor dynamics : 

X= A prj(x- X ) +A ref X + B c r (3.1) 

where A prd is a Hurwitz matrix. 

Using relation (2.15), let 

fe = x-x , 

L „ * (3-2) 

\e = x-x 

represent the tracking and prediction errors, respectively. Subtracting (2.14) from (2.13) and (2.13) from (3.1), 
results in the tracking and prediction error dynamics. 
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(e = A ref e + BA(AK T x x-A& T d (x p )) 


e = A prd e-BA(AK T x x-A® T d ® d 


to) 


(3.3) 


Our design approach is Lyapunov based, and because of that in order to derive stable adaptive laws, consider the 
following Lyapunov function candidate: 

V(e, e, A K x , A® d ) = / P ref e + e T P prd e + trace (AK T x r; 1 AK x a)+ trace (A® T d Tf A0 rf A) (3.4) 

where T x = T t x > 0 and T Qj = r£ > 0 denote constant rates of adaptation, while P ref = > 0 and P prd = P T prd > 0 

are the two unique solutions of the two algebraic Lyapunov equations, correspondingly. 

J Aef ^ref + Kef Aref ~ ~8ref ^ ^ 

| Ard Pprd + P prd A p rd = ~Q prd 

with Q ref = Ql ef > 0 and Q prd = Q T prd > 0 . Time derivative of the Lyapunov function V ( e , e , A K x , A0 rf ), evaluated 
along trajectories of the error dynamics (3.3), yields 

V = Q ref e-e T Q prd e + 2(e T P ref -e T P prd )B a(aK t x x - A0 r , *,(*,)) 


+2 trace ( a/^ T x ' K x A )h 2 trace ( a©^ r 9 ‘ ® d a) 

Define the training error signal: 

e=B T {P ref e-P prd e) 

Applying the well-known trace identity (valid for any two co-dimensional vectors a and b ), 

a T b = trace (ba T ) 

the time derivative in (3.6) becomes: 


V = -e T Q ref e-e T Q ] 


Cprd 1 


+2trace(AATj xe T A)+2trace(AAT x r r; 1 K x A )- 2 trace ( a©^ <S> d (x p )e T a)+ 2trace(\0^ ® d a) 

Collecting similar terms, yields 

v = -e T Q ref e-e T Q prd e 

+2 trace ( aATJ k x + *e r 2 trace (a©" & d - <t> d (x p Y Jv j 


(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 


Consequently, choosing the following PMRAC laws, 


K x = - T x xe T 


(3.11) 


implies that 

V (e, e, A K x , A® d )= -e T Q ref e - e T Q prd ezO (3.12) 

which in turn, proves Uniform Ultimate Boundedness (UUB) [3] of (e, e , A K x , A©^ ). Moreover, both the tracking 
and prediction errors are square integrable: e, e£L 2 . Furthermore, since then x ref , then consequently, 

xEL w . This implies jcEI m . Since ideal parameters (A^,©^) are constant and their estimation errors 
(A K x9 AS d ) are bounded, then their estimated values are bounded as well, that is: K x , ® d EL O0 . Since components 
of the regressor vector (x p ) are piece-wise continuous functions of x p and the latter is bounded, the 

components themselves are bounded. Hence, uEL t n and x, x EL„. Thus e^eEL^ which implies that VEL^. 
Therefore, V is a uniformly continuous function of time. Since V is lower bounded, V 0 , and V is uniformly 
continuous then V tends to a limit, while its derivative V tends to zero. This fact is known as Barbalat’s Lemma, 
[1], [2], [3]. Consequently, the tracking and prediction errors, e and e, tend to zero asymptotically, as t -> oo . 
Moreover, since Lyapunov function (3.4) is radially unbounded, the asymptotic convergence of the errors is global. 
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In other words, the closed-loop tracking and prediction error dynamics are globally asymptotically stable. The 
tracking problem is solved. The end result is summarized below. 

Theorem 3.1 

Consider the uncertain system dynamics in (2.9), operating under PMRAC controller (3.11), with the state 
predictor (3.1). Suppose that the matching relations (2.10) hold. Let the reference model (2.14) be driven by some 
bounded and possibly time-varying reference command r (?). Then: 


a) All signals in the closed-loop system, 

* = Ax + B A (K T X x - @ T d <t> d (x p )+ 0^ 3> d (x p ))+ B c r(t) 

Kef =A ref X ref + B c r(t) 

■ x-A prd (x-x)+A^.x+B c r( t) (3.13) 

K X =~T X x ((* - x ref j P ref -(x-xj P prd y 

= r e, ( X P )(( X - V ) P ref-( X ~ X f P pri ) 

are uniformly bounded in time, where symmetric positive definite matrices P ref and P prd are the two unique 
solutions of the two algebraic Lyapunov equations in (3.5), respectively. 

b) Tracking error e = x - x ref is UUB, square integrable, and globally asymptotically stable, lim||e(f)|| = 0 . 

c) Prediction error e = x - x is UUB, square integrable, and globally asymptotically stable, lim||e(?)|| = 0 . 


Remark 3.1 

The developed PMRAC design methodology can be applied to the original system (2.1), and without adding 
integrated output tracking errors. In this case, the control input is: 

u = K T x x + K?r-Q T d * d (x p ) (3.14) 

where K r E R mxm is the estimated feedforward gain. The system closed-loop dynamics become: 


x = A,x + B A 


(k x -K x )x + (k r -K r Jr-(e„ -0 Je„ (x p ) 


+ BAK ; r 


= A r4 x + BK (hKl x - A0> rf (x p ))+ B c r 


(3.15) 


where K r E:R mxm is the unknown / ideal feedforward gain. The relation B c =BAK T r is the required matching 
condition which guarantees existence of the ideal feedforward gain. The corresponding PMRAC laws are: 

K = -r, * ((* - V I P ref -(*-*/ P P r d ^ 

i r =-r r r^x-x ref ] P ref - (x-xj P prd y (3.16) 

^ = r e, (* P )((^ - x ref J P re/ -(x-xj P prd ) 
where T r = T T r > 0 is the adaptation rate matrix for the feedforward adaptive gain K r . 


Remark 3.2 

Comparing MRAC laws (2.17) with PMRAC laws (3.11), it is evident that the presence of the state predictor (3.1) 
adds the low-pass filtering effects of the prediction dynamics to the direct MRAC laws. 
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Remark 3.3 

PMRAC design extension to MIMO systems with non-parametric uncertainties is straightforward and can be 
accomplished by using well-known in adaptive control robustification methods [1], [2], [4], [5], such as: a) Dead- 
zone, b) a - modification, c) e - modification, and d) Projection Operator. 


IV. Simulation Example 

In this section, we will demonstrate PMRAC benefits using a simplified model which is representative of aircraft 
dynamics. Through simulation evaluations, we compare closed-loop tracking performance of the proposed PMRAC 
and the “classical” MRAC controllers. Our focus is on the transient performance of these two controllers. We show 
that in transient, PMRAC yields less oscillations than MRAC. We also show that PMRAC adaptation rates can be 
increased substantially and without inducing unwanted oscillations that would otherwise be present if using MRAC. 
We immediately note that these claims about transient behavior represent numerical observations only. Efforts to 
obtain formal quantification of transient dynamics provided by these two adaptive controllers constitute an ongoing 
research and will be reported at a later time. 


Our simulation model is chosen to represent longitudinal dynamics of an aerial vehicle, such as an F-16 aircraft. 
Neglecting the effects of gravity and thrust, aircraft longitudinal (short period) dynamics can be written in matrix 
form [6], 


la) 


(Z 

Z\ 



Z, \ 

a_ 

V 

1 + -^- 

V 

la) 

| + 

O 

V 

\ a 

Af 
q / 

K $ J 

x ^ 


\ 6 * / 


^pBL 


(4.1) 


where a is the aircraft angle of attack (AO A), q is the pitch rate, S e is the elevator deflection (control input), V is 
the trimmed (constant) airspeed, (Z a ,Z q ,Z d ) and (M a ,M q9 M s ^ are partial derivatives of the aerodynamic 


vertical force Z and the pitching moment M, with respect to ( a 9 q,d e ) respectively. Numerical values for the 

vehicle aerodynamic derivatives were taken from [6, Example 5.5-3, Table 3.4-3]. These data represent an F-16 
aircraft trimmed at: 


V T = 502 / 1/ sec, Alt=0 ft @300 lb ft GG035~, ca = 2. 11 deg 
The resulting open-loop system matrices are: 


—1.0189 

0.9051 \ 

I, B =1 

1 -0.0022 \ 

v 0.8223 

-1.0774 j 

p 

^-0.1756 / 


where a is in radians, q is in radians / second, and 8 e is in degrees. 


(4.2) 


Three types of matched uncertainties are added to the system: a) linear-in- state uncertainty K T x pert x p , b) control 

effectiveness constant uncertainty A > 0 , and c) nonlinear-in-state uncertainty in the form of (2.2). Numerical 
values for the uncertainties were arbitrarily chosen as: 

( a - a cf 

K l pPat = (-4.6839 -6.1373), A = 0.5, d(x p )= d(a)=0.5e (4.3) 

° 7t 

where the center of the Gaussian was set to a =2 , and its width was o = 0.0233 . This particular selection of 

c 180 

numerical values for K T x pert x p and A is equivalent to 50% increase in the static stability M a , 50% decrease in the 

pitch damping M q , and 50% decrease in the elevator effectiveness M & . These changes imply that the vehicle 

became 50% more statically unstable , lost 50% of its pitch damping ability, and the aircraft controllability decreased 
by 50%. Such drastic and perhaps unrealistic changes were motivated by our intent to demonstrate the effectiveness 
of the proposed PMRAC design methodology. 
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The system total matched uncertainty versus AOA was calculated at q = 0 , and is shown in the figure below. 

Total Matched Uncertainty 



Figure 4.1 : Total Matched Uncertainty versus AOA, at q = 0 


With all the uncertainties included, the aircraft dynamics take the form of (2.1), 

x p -(A pJ , l+ B p AK T Xppen )x p + B p A(u + d(x p )) (4.4) 

V 

where x p = (a qj is the system state vector, A > 0 is the uncertainty in the elevator effectiveness, and the state- 
dependent function d ) represents unknown nonlinear increments in the vehicle aerodynamic Z-force and the 
pitching moment M. The vehicle pitch rate q is selected to represent the system controlled output. 

y - q - (0 l)*, = C p x p (4.5) 

Thus, the control goal is to track any bounded time- varying pitch rate command, 

r = <lcmi (0 ( 4 -6) 

in the presence of the system uncertainties. 


Adding integrated pitch rate tracking error dynamics, 
results in the extended open-loop dynamics, which is in the form of (2.9). 


(4.7) 


V 

a 


(0 
0 

0 M. 
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V V 
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\ 

0 \ 



+ 



V 


M* 

/ 

\ , s ‘) 


a (° *;,.») 


e i\ 

a 

q 


l o \ 


m a 


\ 


d^ + d (a) 

\ » 


/- 1 \ 

0 

1° 

B c 


q.cnui 


(4.8) 


^ = (0 0 l)x 

SH ^ 7 v 

y c 

Conventional MRAC architecture was designed using adaptive laws (2.17). The regressor vector O^(a) 
consisted of 11 a - dependent and evenly spaced Gaussians (RBF-s). RBF centers were placed at [-10:2:10] 
degrees of AOA, and all RBF widths were set to o = 0.0233 . Reference matrix A ref was defined using LQR 
method with 
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(4.9) 


Q = diag ([200 0 0]), R = 1 

The corresponding baseline LQR feedback gains were calculated as: 

k Ibl = (14.1421 3.3331 9.1263) (4.10) 

Using the baseline LQR control feedback, 

(4.11) 

the reference model was formed with: 

A ref =A BL+ BK T xBL (4.12) 

It is straightforward to verify that such a selection of A ref satisfies the matching conditions (2.10). The eigenvalues 

of the reference dynamics, along with their corresponding natural frequencies and damping ratios, are shown in the 
table below. 


Eigenvalue 

Damping 

Natural Frequency 

n/d 

n/d 

rad/sec 

-1.12e+000 + 6.99e-001i 

8.47e-001 

1.32e+000 

-1.12e+000 - 6.99e-001i 

8.47e-001 

1.32e+000 

-1.47e+000 

1.00e+000 

1.47e+000 


Table 4.1 : Reference Dynamics Eigenvalues 


In essence, the reference model dynamics were chosen to represent the baseline closed-loop system operating 
under the baseline LQR controller. Figure 4.2 shows the baseline system response to a series of commanded pitch 
rate doublets, along with the elevator deflection and rate that were required by the LQR controller, while tracking 
the pitch rate doublets, without system uncertainties. As seen from the figure, the elevator deflection and rate are 
well within acceptable limits. 


Pitch Rate Tracking Performance 



Figure 4.2: Baseline Closed-Loop Response to Pitch Rate Doublets without Uncertainties 

With only the baseline controller operating and with uncertainties turned on, closed-loop system tracking 
performance degradation can be clearly observed from the data that are shown in Figure 4.3. However, both the 
control input and its rate remain small. 
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Pitch Rate Tracking Performance 





Next, conventional MRAC laws were constructed by solving algebraic Lyapunov equation (2.18) with 

Qref “ diagflo.l 1 800]) (4.13) 

Rates of adaptation were chosen as 

r x = diag ([l 400 400]), r 0j = 20 (4.14) 

and the total control (elevator deflection) was formed as a sum of the baseline LQR controller and its conventional 
MRAC adaptive augmentation, 

6 e=d'BL+S'AD (4-15) 

where the adaptive component d eAD was defined according to (2.12). We immediately note that since initial 
conditions in the adaptive laws (2.17) can be chosen arbitrarily, the selected adaptive augmentation architecture is 
equivalent to a pure adaptive (no baseline included) controller, whose initial conditions K x (0) are set equal to the 


LQR baseline feedback gains. In summary, the total elevator deflection is given by: 

S e =(K xBL+ K x ]x-@ T d ^ d (a) (4.16) 

With LQR and MRAC controllers operating simultaneously, the reference / baseline closed-loop system tracking 
performance is recovered, with acceptable elevator deflection and rates. The data are shown in Figure 4.4. 


Pitch Rate Tracking Performance 
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Figure 4.4: (LQR + MRAC) Closed-Loop Response in the Presence of Uncertainties 

To demonstrate benefits of PMRAC on the system transient response and control input, PMRAC laws (3.11) were 
constructed using adaptive rates (4.14) and employing identical RBF-s to the conventional MRAC laws. Rather than 
selecting matrix Q prd in (3.5), we use the insight into the physical system present in the selection of P ref and note 

that the predictor dynamics should be faster than the reference dynamics. While still satisfying the positive definite 
requirement of (3.5), we select: 

4W= 3 4*. Pprd = 3P r ef (4.17) 

With the LQR and PMRAC controllers turned on, the reference / baseline closed-loop system tracking performance 
is recovered once again, but with significantly smoother control input than the conventional MRAC design required. 
The response data are shown in Figure 4.5. 


Pitch Rate Tracking Performance 



Figure 4.5: (LQR + PMRAC) Closed-Loop Response in the Presence of Uncertainties 


The (LQR + MRAC) and (LQR + PMRAC) controllers resulted in comparable tracking performance (Figure 4.4 
vs. Figure 4.5), yet the required elevator deflection and its rate were of higher levels and contained unwanted high 
frequency oscillations in the case of the (LQR + MRAC) design. T o better illustrate the advantage of the PMRAC 
controller, we overlaid the (LQR + MRAC) and the (LQR + PMRAC) tracking responses, both in the presence of 
the system uncertainties. Figure 4.6 shows the data. 

Pitch Rate Tracking Performance 



Time, sec 


Figure 4.6: Comparison of (LQR + MRAC) vs. (LQR + PMRAC) Closed-Loop Responses 
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An alternative way of quantifying control input activity is to use the frequency domain, and to plot their 
corresponding FFT coefficients. Figure 4.7 shows significantly reduced control input activity in the PMRAC 
augmentation design. 


Control Input: abs(FFT Coeffi cients) 



The data confirm the fact that that the proposed PMRAC controller provides improved transient response (when 
compared to MRAC) of the closed-loop system, in the presence of uncertainties. A formal proof of this phenomenon 
remains an open problem. 

As discussed earlier, the better transient performance characteristics provided by PMRAC can be attributed to the 
inclusion of the prediction error e into the adaptive laws. The main difference between a conventional MRAC 
augmentation and the proposed predictor-based MRAC augmentation lies in the fact that while the parameters in the 
conventional MRAC are updated based on the tracking error between the plant and the reference model, in the 
PMRAC case the adaptive laws are also designed to minimize the prediction error in the input-output measurements. 
In addition to the tracking error in conventional MRAC laws, the prediction error in PMRAC provides an extra 
information about the system uncertain parameters. This allows the predictor-based adaptive augmentation to retain 
stability and tracking performance, while potentially improving the robustness and transient performance of a 
conventional direct MRAC augmentation. 


V. Conclusion 

In this paper, we presented a state-predictor based direct adaptive tracking design methodology for multi-input 
dynamical systems, with partially known dynamics. Efficiency of the design was demonstrated using short period 
dynamics of an aircraft. Formal proof of the reported PMRAC benefits constitute future research and will be 
reported elsewhere. 
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